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Spitzer Observations of Gamma-Ray Burst Host Galaxies: A 
Unique Window into High Redshift Chemical Evolution and 

Star- formation 

R. Chary 1 , E. Berger 2 ' 3 ' 4 , L. Cowie 5 
ABSTRACT 

We present deep 3.6 /zm observations of three z ~ 5 GRB host galaxies with 
the Spitzer Space Telescope. The host of GRB060510B, at z — 4.942, is detected 
with a flux density of 0.23 ± 0.04 /iJy, corresponding to a rest-frame V— band 
luminosity of 1.3 x 10 10 L , or « 0.15 L*y tZ=3 . We do not detect the hosts of 
GRBs060223A and 060522 and constrain their rest-frame V— band luminosity to 
< 0.1 L* 5 v,z=3- Our observations reveal that z ~ 5 GRB host galaxies are a factor 
of ~ 3 less luminous than the median luminosity of spectroscopically-confirmed 
z ~ 5 galaxies in the Great Observatories Origins Deep Survey (GOODS) and 
the Hubble Ultra Deep Field (UDF). The strong connection between GRBs and 
massive star formation implies that not all star-forming galaxies at these redshifts 
are currently being accounted for in deep surveys and GRBs provide a unique 
way to measure the contribution to the star-formation rate density from galaxies 
at the faint end of the galaxy luminosity function. By correlating the co-moving 
star- formation rate density with co-moving GRB rates at lower redshifts, we 
estimate a lower limit to the star-formation rate density of 0.12 ±0.09 and 0.09 ± 
0.05 M yr" 1 Mpc~ 3 at z ~ 4.5 and z ~ 6, respectively. This is in excellent 
agreement with extinction corrected estimates from Lyman-break galaxy samples. 
Finally, our observations provide initial evidence that the metallicity of star- 
forming galaxies evolve more slowly than the stellar mass density between z ~ 5 
and z ~ 0, probably indicative of the loss of a significant fraction of metals to 
the intergalactic medium, especially in low-mass galaxies. 
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1. Introduction 

Our ability to measure the star-formation rate density (SFRD) at z > 4 relies al- 
most entirely on either narrow-ba nd surveys which detect strong Lya emitting galaxies 



(IHu et al 



2004 



Nagao et al.ll2007D or deep imaging surveys of UV b right Lyman-break 



galaxies jGiavalisco et alibood : bouwens et alil2~006l : IHu fc Oowielhood ). These surveys, by 



virtue of being flux limited, trace the bright end of the galaxy luminosity function down to 
~ 0.04 L^ uv,z=3- The various observ ations have revealed a decline in the SF RD by about a 
factor of 3 between z ~ 3 and z ~ 6 (IBouwens et al.ll2006l ; iBunker et al.ll2004l ). with much of 
this decline being due to the evolution of the bright end of the galaxy luminosity function. 

More than 90% of the estimated SFRD at these redshifts takes place in sub-L* j uv,z=3 
galaxies. In addition, spectroscopic confirmation of high redshift galaxies relies on Lya 
emission, which is easily obscured by dust. There is now increasing evidence for rapid dus t 
production within ~ 1 Gyr of the Big Bang (jChary et al.l 120071 . 120051 ; iMaiolino et al.l 12004 ). 
As a result, quantifying possible dust extinction corrections and measuring the faint-end 
slope of the galaxy luminosity function is essential for minimizing uncertainties in the high 
redshift SFRD. 

Measurements of the metallicity in typical star-forming galaxies at z > 4 is beyond 
the technological capability of the current generation of instrumentation. The relevant 
rest-frame optical emission lines are very weak and are redshifted to the mid-IR. The al- 
ternative approach of studying chemical enrichment through damped Lya absorbers (DLAs) 
detected against backgr ound quasars appears to be limited to z < 5 (IProchaska et al.l 12003 



Songaila fc Cowiell2002l ). and is biased towards tracing the properties of extended halo gas, 



which at lower redshift, significantly underestimates the disk metallicity. As a result, the ap- 



from z = to z 



Tremonti et al. 


2004 




Kobulnickv & Kewlev 


2004: 


Savaglio et al. 


z > 5, where such information should shed light 


on the initial 



stages of mass build-up and metal enrichment in galaxies. 

Long duration GRBs, by virtue of being associated with the deaths of massive stars 
provide a complementary technique for measuring the SFRD and the chemic al enrichment 
histo r y Swift has revolu tionized this study by detecting GRBs out to z ~ 6 (IGehrels et al. 



2004 



Kawai et al.ll2006l ). Prompt spectroscopy of the bright afterglows has now provided 
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a sample of ~ 20 GRBs over a wide r ange of redsh i fts w i th a wealth of metal absorption 
features arising in the host galaxy (e.g I Jensen et al.l l200ll ; ICastro et al.l 120031 ; iHjorth et al. 



20031 : IVreeswiik et all 12004 ; iFvnbo et al 



20061 ). These observations provide a unique win- 



dow into the metallicity a n d gas column densi t y in star forming env i ronments at high 



redshifts (IChen et al. 



2006 




Price et al. 


2007 



Prochaska et al.ll2007al ). Once the afterglows fade away, deep observations of the field can 
also reveal the stellar mass and star-formation rate of the host galaxies, which can then be 
correlated with the inferred metallicities. 

In order to study the host galaxies of high redshift GRBs and take advantage of the 
diagnostics afforded by GRBs, we present Spitzer Space Telescope 3.6 /zm ob servations of 



the hos ts of three GRBs at 4 < z < 5. Building on the constraints provided by lBerger et al. 
( I2006al ) on the host galaxy of GRB 050904 at z = 6.295, we discuss the nature of the host 



galaxies, redshift evolution of the luminosity-metallicity relation and provide an independent 
measure of the high redshift SFRD for comparison with estimates from Lyman-break galaxy 
samples. Throughout this paper, we adopt a Qm = 0.27, Q\ = 0.73, H =71 km s -1 Mpc -1 
cosmology. 



2. Observations 



As part of Spitzer program GO 20000 (PI: Berger) we observed t he fields of GRB s 
060223A (z = 4.406; Berger et al., in prep.), 060510B (z = 4.942; IPrice et all bo07h . 
and 060522 (z = 5.110; Berger et al., in prep.) with the Infrared Array Camera (IRAC; 



Fazio et al.ll2004l ) in the bandpasses centered at 3.6 and 5.8 /im (Table 1). The observations 
were undertaken between September and November 2006, after the afterglows associated 
with the GRBs had faded below the detectability threshold. As shown in Table 1, the GRB 
fields lie in regions with "low" to "medium" -level zodiacal background and cirrus of 13—28 
MJy/sr at 24 /im on the date of the observations. We used 100 s integrations with about 130 
medium scale dithers from the random cycling pattern for total on-source integration times 
of ~ 13000 s at each passband. The nominal 3a point source sensitivity limits are 0.26 and 
2.4 //Jy, respectively. 

Starting with the S14.4.0 pipeline-processed basic calibrated data (BCD) sets we cor- 
rected the individual frames for muxbleed and column pull down using software developed 
for the Great Observatories Origins Deep Survey (GOODS). Due to the presence of bright 
stars in the field, many of the frames at 3.6 jum also showed evidence for "muxstriping" . This 
was removed using an additive correction on a column by colum n basis. The processed BCD 
frames were then mosaiced together using the MOPEX routine (IMakovoz fc Khanll2005l ) and 
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drizzled onto a 0.6" grid. Astrometry was performed with respect to the brightest 2MASS 
stars in the field which showed a peak-to-peak astrometric uncertainty of 0.2" at 3.6 fim. 

The location of the GRB hosts was determined by aligning the Spitzer images against 
images of the afterglow from the Swift UV/optical telescope (060223A and 060522) and the 
Gemini Multi-Object Spectrograph on the Gemini-north 8-m telescope (060510B). For the 
latter, the astrometric uncertainty is 0.09" in each coordinate, while using the UVOT images 
we obtain an astrometric uncertainty of about 0.6". 

All three GRB locations show the presence of nearby (3") brighter galaxies: GRB 060223 A 
has two sources with flux densities of 7.2 and 9.1 /xJy at distances of 1.9" and 2.3" from the 
GRB position; GRB 060510B has a source with a flux of 6.1 /xJy about 3.1" from the GRB po- 
sition; and GRB 060522 has a source with 0.72 /iJy located 1.6" away from the burst position. 
This is not unexpected given the high source density in deep IRAC images. We subtracted 
the contribution of these sources, in order to obtain the strongest possible constraints on the 
flux from the GRB host galaxies. 

Photometry at the position of the host galaxies was performed in fixed circular aper- 
tures of 1.2" radius with appropriate beam size corrections applied as stated in the Spitzer 
Observer's Manual. We clearly detect a galaxy coincident with the position of GRB060510B 
with a flux density of 0.23 ± 0.04 fiJy at 3.6 /mi, and a 3a upper limit of 2.4 /iJy at 5.8 jum 
(Figure 1). For GRBs 060223A and 060522, due to blending from nearby brighter sources 
and the residual effects from muxbleed, we are only able to provide 3a upper limits to the 
flux of the host galaxy (Table 1). 



3. Luminosity and Metallicity of GRB Hosts 



Of the 4 GRB host galaxies at z ~ 5 observed in this program at 3.6 and 5.8 /im, 
(including GRB 050904), only GRB060510B is clearly detected. The observed 3.6 /im flux 
densities/limits for these galaxies correspon d to rest-frame V— band luminosities of ~ 0.15 
L*,v,z=3, where L* ) v, z =3 is about 8 x 10 10 L (IMarchesini et al.l 120071 ; IShapley et al.ll200ll ). It 
is illustrative to compare the properties of GRB hosts with the field galaxy population at 
similar redshifts. 



The GOODS fields have extensive spectros copy of galaxies at high redshift (IVanzella et al 



20051 ; IVanzella et al.ll2006l ; IVanzella et aU2007l ) . There are 275 Lyman-break galaxies in both 
the GOODS fields which are classified as V— band "dropouts" i.e. z ~ 5. The magnitude 
l imit of the GOODS o ptical observations imply that they are brighter than ~ 0.2 L* uv>=3 
(jGiavalisco et al.l 120041 ) . Of these, ~20% have spectroscopic redshifts while ~30% are indi- 
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vidually detected with IRAC. At higher redshifts, z ~ 6, it has been shown that galaxies 
which are individually undetected with IRAC appear to harbor a younge r stellar populat ion 
and have a factor of 10 lower stellar mass than IRAC detected galaxies (jYan et al.l 120061 ). 



As shown in Figure 2, GRB host galaxies are factors of 2 — 3 times fainter than the 
median V— band luminosity of galaxies which have spectroscopic redshifts of 4.5 < z < 5.5 in 
the GOODS field. Furthermore, the luminosities are co mparable to the rest-frame V— band 



lumin osity of GRB hosts studied at lower redshifts (e.g.; IChary et al.ll2002l ; iLe Floc'h et al. 



20031 ). This suggests that GRB host galaxies are unlike the luminous end of the star-forming, 
Lyman-break galaxy population which have had about a factor of 10 increase in their stellar 
mass between z ~ 5 and z ~ 1. They are more typical of the blue, faint end of the galaxy 
V— band luminosity function, a population for which it is difficult to measure redshifts or 
metallicities, in the absence of GRBs, due to their inherent faintness. 

GRB host galaxies at z ~ 0.5 — 3, which have extensive multi-wavelength data, show 
clear evidence for very high specific star formation rates indicating an on-going starburst 
(jCharv et al.l 120021 ; IChristensen et all 12004 ; ICastro Ceron et al.l 120061 ) . We do not yet have 
constraints on the star formation rates in the z ~ 5 host galaxies presented here, due to thei r 
intrinsic faintness in the rest-frame UV (see e.g. iFruchter et al.ll2006l ; iJakobsson et al.ll2005l ). 
However, spectroscopy of the afterglows by iPrice et al.l (120071 ) and Berger et al. (in prep.) 
has revealed a wealth of absorption lines which have been used to derive the metallicity and 
gas column density in the vicinity of the burst. 

Absorption spectroscopy of the three bursts presented here have yielded neutral hydro- 
gen gas densities in their host galaxies of: log[iV(HI)] = 21.6 ±0.1 (060223A), log[iV(HI)] = 
21.3 ± 0.1 (060510B), and log[iV(HI)] = 21.0 ± 0.3 (060522). Thus, all three systems are 
clearly PL As , with c olumn densities near t he median of the distribution for GRB-PLAs 
(IBerger et al.l l2006bl ; IJakobsson et al.l 120061 ). In addition, the metallicities of the GRB 
060223 A and 06051 0B systems have been determined from the detection of weak metal lines. 
For GRB 060522 the signal-to-noise of the spectrum is too low to clearly identify any metal 
lines and an estimate of the metallicity is thus not possible. In the case of GRB 060223A, we 
find an upper limit on the column density of S II of log[iV(SiII)] < 15.3, leading to a metal- 
licity of [S/H] < -1.45. The non-detection of Fe IIA1608 leads to a limit of [Fe/H] < -2.65, 
but we stress that iron can be heavily depleted onto dust grains. From the detection of the 
Si IIA1304 line we find log[iV(SiH)] ~ 15.3, and hence [Si/H] » —1.8. As in the case of iron, 
silicon is also strongly depleted, so we conclude that the metallicity of the GRB060223A 
PLA is in the range of ~ -1.8 to ~ -1.4. For GRB060510B, we use the S IIAA1250, 1253 



lines to measure log 



see also Price et al. 



iV(SiII)1 = 15.6 ± 0.1, and hence a metallicity, [S/H] = -0.85 ± 0.15 



20071 ). 
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The metallicity estimates of the GRB hosts along with their rest-frame B— band lumi- 
nosities (assuming a B — V color of 0, typical of star- forming galaxies) are shown in Figure 
3. Also shown for comparison are the metallicity-luminosity relationships for different sam- 
ples of field galaxies. Despite the one detection and two limits for the luminosity of the 
host galaxies, the figure shows that the redshift evolution of metallicity at a fixed B— band 
luminosity that is seen between < z < 2, clearly extends out to z ~ 5. 

The chemical enrichment of galaxies is directly related to their past history of star- 
formation since supernovae and stellar winds are responsible for recycling the products of 
nucleosynthesis back into the interstellar medium. The stellar mass density is the time 



mass density (p*) (e.g 



Dickinson et al. 


2003; 


Yan et al. 


2006; 


Stark et al. 


2007; 


Charv et al. 



20071 ). with the redshift evolution of the metallicity (Z), we can search for evolution of the 
stellar initial mass function and assess the role of feedback in the build-up of galaxies. The 
Spitzer observations of the hosts are crucial, since they enable metallicity comparisons to be 
made at a fixed rest-frame V— band luminosity, over a wide range of redshifts. 

Due to the fact that we have constraints on the V— band luminosity and metallicity of 
only one z ~ 5 GRB host, we make the assumption that the median metallicity at each 
redshift, is that of a galaxy which has a similar luminosity as the GRB host. This is not 
an unreasonable assumption. Within the observational uncertainties, the slope of the mass- 
metallicity relation appears to be invariant between z ~ and z ~ 2 (Erb et aJj|2006j ). The 
metallicity values are obtained by effectively making a vertical cut at —20.8 mag in Figure 
3 and are determined to be — 0.85±0.15, — 0.35±0.1 and 0.33±0.1 dex at redshifts of 5, 2.3 
and respectively. The average estimated p* at these redshifts are 1.4, 6 and 56 in units of 
10 7 M Mpc~ 3 (See references above). We performed a Monte-Carlo analysis to obtain the 
best fit between Z(z) and p*(z). 

Star-forming galaxies which fall on the local mass-metallicity relat ionship, show a scatte r 
of ~0.1 dex at bright luminosities and ~0.2 dex at faint luminosities (ITremonti et al.l 120041 ) . 
We use a random number generator to offset the stellar m ass density and metal licity by the 
observed scatter from the mean values quoted above (See iDickinson et al.ll2003l . Table 3 for 
the range in stellar mass density). We fit for the relation between Z(z) — p*(z) and repeat 
the process 10000 times. We find that dZ/dp* appears to be invariant between < z < 5 



and that Z(z) oc p*(z) 



0.69±0.17 



This suggests that the chemical enrichment of star-forming 



galaxies takes place at a slower rate than the build up of stellar mass. This is presumably 
due to the loss of metals from low-mass galaxies by outflows and stellar winds, an effect 
which is primarily re sponsible for the mass-metallicity relation seen i n the local Universe 
(ITremonti et al.ll2004l ) and z ~ 2 Lyman-break galaxies (Erb et al.ll2006l ). However, alternate 
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mechanisms such as depletion of metals onto dust grains cannot be ruled out at this time. 



There is the possibility of a selection effect in this analysis. If long duration GRBs 
arise in collapsars, they might preferentially be in low-metallicity galaxies. As a result, it is 
possible that GRB hosts have a lower metallicity than the average field galaxy of the same 
rest-frame optical luminosity. Although GRB hosts appear to be have low luminosities in 
the rest-frame UV and V— band, the observational evidence does not indicate that the hosts 
have an unusually low-metallicity for th eir luminosity M e tallicity of GRB host galaxies 
appear to span the range 0.1—1 Z so i ar (e.g lBerger et al.ll2007l . l2006bl ; iProchaska et al.ll2007bl ) 
and some of the hosts have even b een found to be associated with dusty, infrared luminous 
galaxies (e.g 



Le Floc'h et all 120061 ) 



Nevertheless, we assess the reliability of our derived Z(z) — p*(z) relation by considering 
a bias in the metallicity of GRB environments. If we assume that the metallicity of the 
GRB environment is higher by >0.3 dex compared to the mean metallicity of a galaxy at its 
luminosity, it implies that the mean metallicity at z ~ 5 for a field galaxy at the luminosity 
of the GRB host is -1.15±0.15. The best fit relation to the three points is then consistent 
with an exponent of unity i.e. Z(z) oc p*(^) ' 85±0 ' 19 but has a worse x 2 . The corollary is 
that if GRB hosts were biased by 0.3 dex towards lower metallicities, compared to the mean 
metallicity of a galaxy at its luminosity, the best fit relation is Z{z) oc p*(z) a52±a16 which is a 
larger deviation from unity. Furthermore, if there were a bias in GRB host metallicities, the 
slope of the Z(z) — p*(z) relation derived above at a fixed B— band luminosity, would have a 
different value between 2 < z < 5 and < z < 2 due to the fact that the < z < 2 relation 
is determined from star-forming galaxies while the 2 < z < 5 relation is derived from GRB 
hosts and Lyman-break galaxies. This is inconsistent with our fits, although larger samples 
of GRB hosts are needed to eliminate suggestions of bias. 

Detection of individual GRB hosts at high redshifts is likely to remain difficult, due to 
their intrinsic faintness. There is a clear need for homogeneous infrared surveys of GRB host 
galaxies which will enable stacking to be performed as a function of metallicity, gas density 
and rest-fr ame ultraviolet prop erties. Within our sample, GRB050904 is a marginal IRAC 
detection (IBerger et al.ll2006al ). while GRB060223A is dominated by detector systematics. 
As a result, we are unable to provide additional constraints using stacking. Observations 
of a larger sample of GRB hosts, such as those currently being targeted in Spitzer program 
GO4-40599 (PI: Chary), will allow the luminosity- metallicity relation to be measured at high 
redshift and lead to a better understanding of the faint end of the galaxy luminosity function, 
a regime which is currently inaccessible even through ultradeep surveys like GOODS and 
the UDF. 
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4. Evolution of the Star-formation Rate Density 



It is now well known from various mid- infrared, far-infrared and submillimeter surveys 
that the star-formation rate density at z ~ 0.5 — 2.5 is dominated by inf rared luminous 

and L IR /L UV ~ 10 — 100 (e.g iTakeuchi et al. 



2005; 


Burearella et al. 


2006: 


Charv & Elbaz 


2001 



At z > 3, current long- wavelength sur- 
veys, due to their limited sensitivity, are unable to detect galaxies which harbor the bulk of 
the star-formation. Thus, rest-frame ultraviolet observations of galaxies are the only avenue 
for probing star-formation at high redshifts. 

The primary uncertainties associated with quantifying the SFRD at 2 > 3, are the 
contribution from galaxies at the faint end of the UV-luminosity function and dust correc- 
tions. Since sub-L* j uv,z=3 galaxies contribute ~90% of the SFRD, measurement of the faint- 
end slope of the ultraviolet luminosity function, where completeness correction s and surface 



brightness dimming issues are significant, needs to be undertaken carefully (jSteidel et al. 



19991 ; iBouwens et al.ll2006l ). Similarly, if extinction were a significant issue, the galaxies that 
dominate the star-formation rate density would be UV-faint or undetected in magnitude 
limited rest-frame ultraviolet surveys. GRBs are relatively insensitive to these limitations. 
If the GRB rate density were correlated with the co-moving star-formation rate density at 
lower redshifts, where cross- calibration between the UV and IR are in broad agreement, 
measurement of the GRB rate densi ty at z > 3 coul d provide an independent pathway to 
quantifying the SFRD (see also e.g.; iPrice et al.ll2006h . 



The three parameters which are most likely to dominate the calibration between GRBs 
and the star- formation rate density are the evolution of metallicity with redshift, evolution of 
the initial mass function of stars and identification and spectroscopic follow up of the GRB 
afterglow. If long duration GRBs were to preferentially occur in low metallicity environments, 
the increase in the average metallicity of the Universe with decreasing redshift would result 
in a higher SFR/GRB-rate ratio at low redshift. Similarly, evolution of the stellar initial 
mass function from a "top-heavy" to a Salpeter mass function with decreasing redshift 
would increase the SFR/GRB-rate ratio at low redshift. On the other hand, the detection 
efficiency and spectroscopic completeness of GRBs should be increasing with decreasing 
redshift, implying a lower SFR/GRB-rate ratio at low redshift. 

Calibrating each of these parameters individually is challenging at the present time, 
partly because the relationship between GRB rate and environment is not well known and 
due to the fact that observational selection effects cannot be quantified. Therefore, we need 
to rely on empirical comparisons between known star-formation rate estimates and GRB 
rate densities to assess GRBs as a star-formation rate indicator. This empirical comparison 
can be optimally done at z < 3 since in this redshift range the star-formation rate, including 
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the dust obscured component, has been accurately determined from deep mid- infrared and 
submillimeter surveys. 



We use the star-formation rate density at z < 3 from IChary &: Elbaa (120011 ). We 
distribute the 52 Swift GRBs with spectroscopic redshifts into redshift bins and divide by 
the co-moving volume in each redshift bin. We also correct for the time dilation to estimate 
the comoving GRB rate density over the ~ 2 year Swift lifetime. The redshift bin at z < 0.5 
is omitted since the GRB rate density appears to be anomalously high compared to the 
rapidly evolving star-formation rate density. We find that within the uncertainties, the rate 
density of GRBs with spectroscopic redshifts in the range 0.5 < z < 3 is constant at a value of 
(3.7±1.1) x 10 -11 Mpc" 3 yr" 1 . This can be compared with the extinction-corrected comoving 
star-formation rate density in the same redshift range which is in the range 0.12—0. 25 M 
yr" 1 Mpc -3 and has an average value of ~ 0.2 M Q yr" 1 Mpc" 3 ( IChary fc Elbaa 1200 ll ). 



Since these two independent rate densities are relatively constant in the 0.5 < z < 3 
range, we can tentatively make the assumption that the SFR/GRB-rate is constant (Figure 
4). The ratio of these two rates implies: 

SFRD = GRB rate x (5.2 ± 2.3) x 10 9 , (1) 

where SFRD is the extinction-corrected star-formation rate density in M yr" 1 Mpc" 3 and 
GRB rate is in units of Mpc" 3 yr" 1 . 



Using our derived calibration, and the measured GRB rate densities at 4 < z < 5 
and 5 < z < 7 of (2.4 ± 1.2) x 10" 11 and (1.8 ± 0.9) x 10" 11 Mpc" 3 yr" 1 , respectively, 
we infer a net star- formation rate, corrected for extinction, of 0.12 ± 0.06 and 0.09 ± 0.05 
M yr" 1 Mpc" 3 at z ~ 4.5 and z ~ 6 respe ctively. These estimates are systematically 
higher than those derived by iPrice et al.l (120061 ) by factors of 3 — 5. The iPrice et al.l (120061 ) 



estimates were calibrated at z ~ 3 where neither the completeness correction factor for the 
faint end of the UV luminosity function nor the dust extinction correction are reliably known 
while deep Spitzer mid-infrared surveys have confirmed the dominant c ontribution of infrared 



luminous galaxies to the star-formation rate density at 0.5 < z < 3 (ITakeuchi et al.l 12005 



Daddi et al.l 120071 ). The fact that the GRB rate density is almost flat between 0.5 < z < 6, 
while parameters such as the detection efficiency and spectroscopic completeness should be 
decreasing with increasing redshift, implies that the measured GRB rate density provides at 
least a lower limit to the star-formation rate density. 

It is illustrative to compare this star formation rate est imate with those from deep 
rest-frame ultraviolet surveys at z > 4. iGiavalisco et al.l (120041 ) derive a star-formation rate 
density at z ~ 4 of 0.02 M yr" 1 Mpc" 3 when integrating to 0.2 L*,uv,z=3- After application 
of an extinction correction of Ay = 0.45 mag, based on the extinction properties in local 
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starburst galaxies, they esti mate the total s t ar-for mation rate density at z ~ 4 to be 0.15 
M yr _1 Mpc -3 . Similarly, iBouwens et al.l (120061 ) derive a star- formation rate density at 



z ~ 6 by integrating the luminosity function of Lyman-break galaxies in the UDF and other 
deep fields. Integrating the UV luminosity function down to 0.2 £*,uv,z=3 results in a value 
of 1.3 x 10~ 2 M yr" 1 Mpc~ 3 while the integral to 10 7 L© yields a SFRD of 0.04 M yr" 1 
Mpc~ 3 . Application of an extinction correction, inferred to be about A\jy=0A5 mag at 
z ~ 6, to this latter number, implies a SFRD of 0.06 M Q yr _1 Mpc -3 . 

The agreement between the SFRD values estimated from ultraviolet surveys and the 
GRB rate density is reassuring, considering that there have been only 8 GRBs that have 
been spectroscopically confirmed to be at z > 4 (Figure 4). However, the SFRD from 
GRBs primarily traces the faint end of the galaxy luminosity function while the surveys 
are measuring the contribution from the bright end. As a result, a more reasonable SFRD 
estimate requires adding the SFRD contribution estimated from the faint end of the galaxy 
luminosity function, from GRBs, to that from bright LBGs. 

GRB hosts are fainter than 0.2 L*,v,z=3- Based on the UV to V— band flux ratios of 
star-forming galaxies at z ~ 3, it implies that GRB hosts must be fainter than 0.2L* j uv,z=3 
Adding the SFRD from L > 0.2 L*,uv,z=3 galaxies to that inferred from the GRB rate 
density results in an extinction corrected SFRD of 0.27 ± 0.13 and 0.11 ± 0.05 M yr _1 
Mpc~ 3 at z ~ 4.5 and z ~ 6, respectively. If confirmed through a larger statistical sample, 
this is a substantial upward revision suggesting that L < 0.2 L* ! uv,z=3 galaxies contribute 
at least four times as much to the star-formation rate density at z ~ 6 as the bright end 
(L > 0.2 -L*,uv, z =3) of the UV luminosity function. Indirectly, this implies that the faint end 
slope of the UV luminosit y function at z ~ 6 m ust be ~ —1.9, compared to the value of 



-1.73 that was derived by lBouwens et al.l (120061 ). 



GRBs are a powerful tool for measuring the high redshift star-formation rate density. 
In particular, deep Spitzer observations of GRB hosts can reveal the contribution to the 
star-formation rate density from the faint end of the galaxy luminosity function, a regime 
which is inaccessible to deep, rest-frame ultraviolet /near-infrared surveys. Increasing the 
sample of high redshift GRBs will reduce the uncertainties in the star-formation rate density 
unaffected by extinction and through stacking analysis on the host galaxies will help estimate 
the contribution to the stellar mass density from sub-L* galaxies. Comparison between star- 
formation rate estimates from GRBs with those from deep UV surveys will provide better 
constraints on the evolution of dust extinction at high redshift and provide tremendous 
insights into the chemical enrichment of the early Universe. 
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Fig. 1 — Spitzer image of the host galaxy of GRB060510B at 15 h 56 m 29.607 s ,+78°34'12.42" 
(J2000). Image is 12" on a side, North is up, East to the left. The left panel shows the 
processed mosaic while the right panel shows the image with the foreground galaxy 3.1" to 
the East subtracted. 
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Fig. 2. — (Left) Histogram showing the distribution of observed 3.6 /zm magnitudes for 
galaxies in the GOODS fields with spectroscopic redshifts 4.5 < z < 5.5. The solid symbols 
show the brightness of the GRB host galaxies observed in this paper relative to the field 
galaxies. (Right) GRB hosts have rest-frame V— band luminosities which are a factor of ~ 2 — 
3 fainter than field galaxies at similar redshifts and provide a complementary way to study the 
faint end luminosity function of star-forming galaxies. Also shown as th e shaded histogram 
is the V— band lumin osities of GRB hosts at a median redshift of ~1 (I Chary et all I2002 



Le Floc'h et all I2003I ) which indicate that GRB hosts span similar V— band luminosities, 



regardless of redshift. 
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Fig. 3. — Luminosity-metallicity relationship for star-forming galaxies at z < 2 compared 
with the host galaxies of z > 4 gamma-ray bursts. The galaxy data are from GDDS and 
CFRS at z ~ 0.4 - 1 (cir cles; Isavaglio et aD l2005h . TKRS at z ~ 0.3-1 (diam onds; 
Kobulnickv fc Kewlevl 120041) DEEP2 at z ~ 1 - 1.5 (squares; Ishaplev et all l2005h and 



LBGs at z ~ 2.3 (error bars; lErb et al.ll2006l ). The gray lines represent the relation derived 
for z ~ 0.1 galaxies from the Sloan Digitized Sky Survey (ITremonti et al.ll2004l ). GRBs 
provide a unique window into the evolution of the mass-metallicity relation at high redshift 
and indicate that the chemical enrichment of galaxies with redshift occurs at a lower rate 
than the build up of stellar mass, presumably due to the expulsion of metals in low-mass 
galaxies by outflows. 
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Fig. 4. — Star-formation rate density inferred from spectroscopically confirmed long-duration 
Swift GRBs (solid squares). The solid black line is the extinction-corrected star- formation 
rate density inferred at z < 4 from a variety of multiwavelengt h surveys in the mid - 
infrared and submillimeter, which are used to calibrate the GRBs ( jChary fc Elbaa l200ll ). 
The lower hatched region is the extinction uncorrected SFRD fr om rest-frame UV su rveys 
(> 0.04 L^uv^s), in cluding estimates by ISteidel et all (Il999l ) , lYoshida et all (120061 ) and 



Bouwens et all (120061 ). The upper hatched rectangle are these y alues corrected upward for 



reddening using the UV-slope technique by lBouwens et all (120061 ). The SFRD inferred from 
GRBs at z > 4 is consistent, within the significant errors, to the extinction corrected SFRD. 
Follow-up of a larger number of high-redshift GRBs are required to confirm if the higher rate 
density derived from the GRBs is statistically significant. 



Table 1. Spitzer Observations of z ~ 5 GRB Host Galaxies 



GRB 


RA, DEC (J2000) 


Rcdshift 


Date of Observation 


Sky Background a 


Exposure Time 


Flux Density 


in ^Jy 










MJy/sr 




3.6/rni 


5.8/im 


GRB 060223A 


03 h 40 m 49.561 s , -17°07' 48. 36" 


4.406 


2006 Sep 23 


24.8 


130 x100s 


<0.3 


<2.4 


GRB 0605 10B 


15 h 56 m 29.607 s ,+78°34'12.42" 


4.942 


2006 Oct 26 


13.4 


130 x100s 


0.23±0.04 


<2.1 


GRB 060522 


21 h 31 m 44.800 s ,+02°53'10.35" 


5.110 


2006 Nov 23 


28.4 


138 x100s 


< 0.2 


<2.4 



a Background at 24//m, dominated by the zodiacal light. 



